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Abstract. We have run high resolution N-body models simulating the encounter of a dwarf galaxy with a bright elliptical 
galaxy. The dwarf absorbs orbital angular momentum and shows counter-rotating features in the external regions of the galaxy. 
To explain the core-envelope kinematic decoupling observed in some dwarf galaxies in high-density environments requires 
nearly head-on collisions and very little dark matter bound to the dwarf. These kinematic structures appear under rather re- 
strictive conditions. As a consequence, in a cluster like Virgo ~ 1% of dwarf galaxies may present counter-rotation formed by 
harassment. 

Key words. Galaxies: dwarf- galaxies: interactions- galaxies: kinematics and dynamics- galaxies: structure- methods: N-body 
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1. Introduction 

Dwarf galaxies, the most numerous galaxy type in the Universe 
(Binggeli, Sandage, & Tammann 1988; Ferguson & Binggeli 
1994; Gallagher & Wyse 1994), come in two groups: dwarf el- 
lipticals (dE) and dwarf irregulars (dl). These low mass systems 
have similar stellar distributions, showing exponential surface 
brightness profiles with similar central surface brightness and 
scale lengths (Lin & Faber 1983). They also follow the same 
luminosity-metallicity relation (Skillman et al. 1989; Richer & 
McCall 1995). But, dEs and dls differ in the gas content, dls 
being gas-rich galaxies while dEs are devoid of warm gas. 

Whether dls and dEs are evolutionary linked, and what such 
links may be, is a matter of current debate. Both internal and ex- 
ternal processes have been proposed. Among the internal pro- 
cesses, kinetic energy from supernova explosions can sweep 
the gas and turn dl galaxies into dEs (e.g. Dekel & Silk 1986; 
De Young & Gallagher 1990). However, such mechanism can- 
not explain recent evidence that the velocity gradients in most 
dE galaxies in Virgo cluster are generally different from those 
of dls (van Zee, Skillman & Haynes 2004). 

Environmental processes may be suspected given that dE 
galaxies are the most strongly clustered galaxy type (Ferguson 
& Sandage 1989), hence encounters with other galaxies could 
be important in their evolution (Aguerri et al. 2004, 2005 and 
references therein), van Zee, Skillman & Haynes (2004) pro- 
pose interaction-driven gas striping as the main mechanism for 
the transformation of dls into dE galaxies. Dwarf galaxies may 
also suffer transformations due to fast interactions with large 
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cluster members and to their motion in the steep cluster poten- 
tial, a process known as harassment (Moore et al. 1996; Mayer 
etal.2001). 

Understanding the fossil records of evolutionary transfor- 
mations in a high-density environment may help sorting out 
the relative importance of the above mentioned processes. In a 
recent paper, de Rijcke et al. (2004) argue that fast interactions 
with giant galaxies may leave a signature in the rotation curve 
of dwarf galaxies, including core-envelope counterrotation. In 
their deep spectroscopy of 15 dE galaxies, they report on two 
cases of dEs showing complex kinematic profiles. FS373, a nu- 
cleated dE(2,N) in the NGC 3258 group, shows a 0.3 kpc (-0.2 
R e ) kinematically-decoupled core (KDC), in corrotation with 
the main body, plus an outer velocity reversal at 2.2 kpc (~1.5 
R e ). FS76, a dEO dwarf in the NGC 5044 group, shows a nu- 
clear corrotating KDC in the inner ~0.2 kpc (0.26 R e ) and an 
asymmetric rotation profile outside 1 R e . 

On the basis of the large velocity dispersions of both 
groups, de Rijcke et al. (2004) argue that mergers are unlikely 
mechanisms for the origin of the complex kinematic profiles of 
the two dEs. Instead, they propose that the velocity reversals 
could result from fast interactions with the dominant galaxies 
in their groups. Their analytic estimates based on the impulse 
approximation show that sufficient angular momentum may be 
exchanged during a fast encounter to reverse the rotation of the 
dwarfs. 

Kinematic decoupled cores (KDC) in bright galaxies have 
been well studied in the past two decades. Most well known 
are counterrotating cores in elliptical galaxies (Bender 1990), 
which have been explained with mergers of unequal ellipticals 
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(Balcells & Quinn 1990), or with mergers of spiral galaxies, 
either with (Hernquist & Barnes 1991) or without (Balcells 
& Gonzalez 1998) dissipative gas dynamics. Some kinematic 
reversals in elliptical nuclei might be the result of projection 
effects on triaxial ellipticals (Oosterloo, Balcells, & Carter 
1994). Counterrotating bulges may originate in mergers of 
dwarf galaxies onto large disk galaxies (Aguerri, Balcells, & 
Peletier 2001). We note that counterrotation may appear in pri- 
mordial collapses as well (Harsoula & Voglis 1998). 

For dwarf galaxies, while mergers can in principle yield 
counterrotation, we concur with de Rijcke et al. (2004) that 
high group velocity dispersions make those highly improbable. 
Other processes, studied for the cases of counterrotation in spi- 
ral galaxies, may apply to dwarfs. These include: scattering of 
stars by a bar (Evans & Collett 1994); the evolution of a polar 
ring in a triaxial halo (Tremaine & Yu 2000); and even simple 
projection effects on triaxial or non-planar stellar disks. 

The harassment hypothesis put forward by de Rijcke et al. 
(2004) offers yet another mechanism for the origin of KDCs in 
dwarf galaxies. In this paper we simulate, using high-resolution 
numerical models, the impulsive encounter of a dwarf galaxy 
with a giant elliptical, in order to perform direct measures of the 
true effect of the encounter on the rotation profile of the dwarf. 
Our aim is to test the validity of the analytical predictions of 
angular momentum exchange during impulsive encounters. A 
similar study was published by Mayer et al. (2001), who mod- 
eled interactions of dls with the halo of the Milky Way to ex- 
plore dwarf evolution by tidal stirring. They do not report on 
any KDC in their models, which include retrograde interac- 
tions. The two studies differ in that Mayer et al. set the satel- 
lites on bound orbits and model the larger galaxy as an exter- 
nal potential, thus neglecting dynamical friction. Our models 
describe fast, hyperbolic, strongly penetrating encounters with 
a live primary galaxy. Our stronger interactions provide more 
stringent tests on the formation of counterrotating cores in fly- 
by interactions. 

In our paper, we first review analytical predictions (§ 0, 
then describe the models (§EJ. Our results (§0} will show that, 
while analytical predictions of angular momentum exchange 
are accurate, the effects on the rotation curves are small: the 
dwarf material that absorbs most of the angular momentum ab- 
sorbs energy as well, and escapes from the dwarf. Envelope 
counter-spinning is attained for nearly head-on collisions only, 
and only if the dwarf does not have a dark matter halo. We 
conclude that fly-by interactions leading to counterrotation at 
R < 1 .5R e may have occurred in about 1 % of the dwarf galax- 
ies in a Virgo-type cluster. 

2. Theoretical expectations 

We here review analytical predictions for the exchange of an- 
gular momentum in a fly-by encounter of a dwarf galaxy with 
a giant galaxy. Our goal is to see under what conditions the 
rotation of a dwarf envelope can reverse sign as a result of 
the fly-by. de Rijcke et al. (2004) made one such prediction 
based on the tidal transfer calculation of Som Sunder, Kochhar 
& Alladin (1990), yielding a global value for the change of the 
rotation of the dwarf in terms of its inertia tensor. We slightly 
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Fig. 1. Velocity change according to the impulse approximation 
(eqs. Q and (0) for a dwarf galaxy flying through a normal 
galaxy with mass ratio 10: 1 on a parabolic encounter with peri- 
center distance close to the half-mass radius of the large system 
(model CeP in Table[l). Solid light curve is the initial velocity 
curve, the dotted-dashed light curve is the velocity change ac- 
cording to eq. [2 the dotted-dashed heavy curve is the velocity 
change according to eq.|2 and the light dashed curve and solid 
dashed curve are the final velocity curves accordingly. 

modify that calculation in order to estimate the velocity change 
as a function of radius in the dwarf, so that diagnostics on coun- 
terrotation can be inferred. 

We use the impulse approximation: we assume that the 
galaxies have no time to reorganize during the encounter, hence 
the energy injection is all kinetic; reorganization occurs after 
the encounter has finished. We further make use of formulae 
based on the tidal approximation (Binney & Tremaine 1987, 
§ 7.2(d)), in which the tidal field has been expanded to first or- 
der. Strictly speaking, the impulse approximation applies only 
if the encounter time is much shorter than the internal cross- 
ing time, and the tidal approximation applies when the impact 
parameter is much larger than a typical radii of the galaxies. 
While such conditions do not strictly apply to some of the mod- 
els we simulate, these approximations provide useful estimates 
of the energy and angular momentum exchange during fast en- 
counters (e.g., Aguilar & White 1985). 

In the tidal approximation, the velocity change A Vi in the 
stars of dwarf galaxy after a fast interaction with a perturber of 
mass Mi scales with the dwarf galactocentric radius R as 

lAVaOOh^J (1) 

where b the impact parameter of the encounter, V the relative 
velocity of the galaxies and G is the gravitational constant. We 
have omitted terms of order 1 that depend on the dwarf's axis 
ratios. A V2(R) linearly increases with R, hence the perturbation 
increases in the outer parts. 

In the case when the spin of the dwarf galaxy is opposite 
to the orbital angular momentum, the change in rotation veloc- 
ity given by eqn.n can lead to core-envelope counterspinning. 
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Fig. 2. Left: Binding energy for the particles of the dwarf galaxy in model CaHH at the end of the fly-by interactions. All particles 
from the dwarf are plotted. Middle: Energy change according to the impulse approximation (eq.[3J- The dashed curve is the initial 
binding energy for particles on circular orbits. The dotted line is the expected energy gained according to eq.[3] and the solid line 
is the final binding energy. Right: Energy change according to eq.0] thin lines provide the correction due to the extended nature 
of the giant, thick lines give the full approximation. Dashed, dotted and solid lines are as before. 



As an application of eqn.^ we show in Figure[2the change in 
the rotation curve expected in a dwarf galaxy as a function of 
the galactocentric dwarf distance scaled with the effective ra- 
dius of the dwarf R e 2- The change after a retrograde, parabolic 
encounter with a 10 times more massive galaxy is shown. The 
initial rotation curve and effective radius of the dwarf are those 
used in the iV-body models, and the pericenter distance is equal 
to the effective radius of the massive galaxy (see §[3}. We find 
that a velocity change as given by eqn. [pleads to a counterro- 
tation radius R a * 1 .5/? e 2- 

The tidal approximation used in Equation [2 considers the 
large system as a point mass. We can generalize eqn.nto take 
into account the extended nature of the large system. Following 
Gnedin, Hernquist & Ostriker (1999) we have: 



GM\ Rneri t 



(2) 



where R pe ri is the pericenter distance and R max is the cut-off 
imposed in our initial model for the large galaxy (see §[3J. The 
curves predicted by this new formula can be seen in Figure ^ 
as thick lines. Now the velocity change is smaller and counter- 
rotation is found at larger radii. 

Equations Hand 13 provide AV2 under the assumption that 
all the transferred angular momentum is kept within the dwarf. 
But stars in the dwarf absorb energy as well as angular momen- 
tum. We can estimate the variation of internal energy per unit 
mass in the secondary at radius R by A.E2 ~ 0.5(AV2) 2 , given 
that, in the impulse approximation, all the energy absorbed is 
kinetic: 



AE 2 (R) 



G M\ 
b 4 V 2 



AE:(*>~^tf[(jH 2 <l -a/0-"] 



(4) 



R 1 



(3) 



tfV 2 " ^Rmax 

where equation [3] comes from Equation [2 and eq. 0] includes 
the correction due to the extended nature of our giant galaxy 
and the term due to the adiabatic corrections, i.e. the correc- 
tions that we have to take into account due to the fact that the 
stars will move during the interaction. Here u is the orbital 
frequency of stars in the satellite and r the length of the tidal 
shock. The exponent -2.5 in fact changes accordingly to the du- 
ration of the tidal impulse relative to the dynamical time of the 
satellite from 2.5 - 3 to 1.5 (see Gnedin et al. 1999) 

Particles that become unbound as a result of the energy in- 
jection carry out part of the angular momentum delivered by 
the interaction. The escape of particles is confirmed by Figure|3] 
which shows the energy budget for one of our models and the 
expressions given in eq.[3]and0] Figure|2(left), shows the dis- 
tribution of binding energy in the dwarf after the interaction, 
for the A^-body realization of model CaHH. A tail of unbound 
particles is seen at R/R e > 3 - 4 in this snapshot which corre- 
spond to the last time step of our model, when the system has 
evolved after the encounter. Middle panel shows the effect of 
the pure tidal approximation (eq.[3}- The figure plots the bind- 
ing energy per unit mass as a function of radius before the inter- 
action (circular orbits; dashed line); the energy change AE2(R) 
from eqn. [3] (dotted line); and the expected final distribution 
of binding energy (solid line). Many particles outside 2 R e be- 
come unbound and escape the dwarf. Failing to account for the 
angular momentum carried away by these particles should lead 
to a large overestimate of the angular momentum delivered to 
the dwarf, and to a corresponding overestimate of the effects on 
the dwarf rotation curve. Figure|2]right panel shows the results 
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Table 1. Initial configurations for fly-by models . 



Mod. 




(6*2, <fe) 


n 


e 


''peri 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


CaH 


10:1 


(180,0) 


20 


1.07 


0.05 


6.77 


CaHH 


10:1 


(180,0) 


20 


1.003 


0.05 


6.64 


CbH 


10:1 


(180,0) 


20 


1.07 


0.10 


4.77 


CbHH 


10:1 


(180,0) 


20 


1.013 


0.10 


4.70 


CbP 


10:1 


(180,0) 


20 


1.00 


0.10 


4.69 


CcP 


10:1 


(180,0) 


20 


1.00 


0.20 


3.32 


CcH 


10:1 


(180,0) 


20 


1.07 


0.20 


3.37 


CdP 


10:1 


(180,0) 


20 


1.00 


0.41 


2.31 


CdH 


10:1 


(180,0) 


20 


1.07 


0.41 


2.35 


CeP 


10:1 


(180,0) 


20 


1.00 


0.83 


1.63 


CeH 


10:1 


(180,0) 


20 


1.07 


0.83 


1.66 


CfP 


10:1 


(180,0) 


20 


1.00 


5.00 


0.66 


CfH 


10:1 


(180,0) 


20 


1.07 


5.00 


0.68 


CeP G 


10:1 


(180,0) 


20 


1.00 


0.83 


1.63 


CeP c 


10:1 


(180,0) 


20 


1.00 


0.83 


1.63 


CeP H 


10:1 


(180,0) 


20 


1.00 


0.83 


1.63 


CeP R 


10:1 


(180,0) 


20 


1.00 


0.83 


1.63 


CcP R 


10:1 


(180,0) 


20 


1.00 


0.2 


3.32 


CbP R 


10:1 


(180,0) 


20 


1.00 


0.1 


4.69 


ReP 


10:1 


(0,0) 


20 


1.00 


0.83 


1.63 


CReP 


10:1 


(135,0) 


20 


1.00 


0.83 


1.63 


CePl 


5:1 


(180,0) 


20 


1.00 


0.83 


1.71 


CeP3 


20:1 


(180,0) 


20 


1.00 


0.83 


1.60 


CePA 


50:1 


(180,0) 


20 


1.00 


0.83 


1.57 



for the extended nature of the primary (light curves) and for the 
full eq.|4](heavy curves) which are in closer agreement with the 
results from our experiments, as we discuss later. 

TV-body fly-by experiments are the ideal tools to determine 
the true effects of the interactions on the rotation curves of 
dwarf galaxies. We address those in the subsequent sections. 

3. The Models 

Our numerical models simulate the fly-by of a dwarf galaxy 
about a massive elliptical. A King (1966) model was adopted 
for the dwarf. The King model has an exponential surface den- 
sity profile over a large radial range, making it a fair approxi- 
mation for dE galaxies (Barazza, Binggeli, & Prugniel 2001). 
The model was built using the GalactlCS code (Kuijken & 
Dubinski 1995). The central density was 14.45 and the cen- 
tral velocity dispersion 0.714 (in internal units, see below). The 
model was then scaled down, using M ~ r 2 (Fish 1966), to 
mass 0.1, and was allowed to relax in isolation for fifty half- 
mass crossing times. During the initial 8-10 crossing times 
the system expanded, probably due to the softening (about 1/5 
of the half-mass radius), then remained in equilibrium for the 
rest of the run. We used the final relaxed model after the 50 
crossing-times as input for our interaction experiments. After 
scaling and relaxation, the central density and velocity disper- 
sion are 10.40 and 0.302, respectively. The effective radius is 
R e 2 - 0.178, while the King concentration index is C = 0.72, 
and the radius that encloses the 99% of the mass is 2.76. We 
imposed a small rotation by reversing the velocities of 80 % of 
the particles chosen randomly, yielding V ma x/co = 0.49. 



In order to compare our initial dwarf galaxy model with 
observed dwarf galaxies we have used the sample by Geha, 
Guhathakurta & van der Marel (2003). We use the rotating dE 
VCC 1947 from their sample as our fiducial initial system. This 
dE has an absolute magnitude My = -16.32 and an effective 
radius of 0.62kpc. According to this number the mean density 
inside the effective radius would be ~ 2M Q /pc 3 (assuming a 
mean value of M/L = 18M Q /L derived from dwarfs in the 
Local Group (Mateo 1998)). Using the units given below, our 
initial system yield a mean density inside the effective radius 
1.5M / pc 3 . Figure[3]top panel compares the rotational support, 
expressed as V(R)/cr(R), for our initial model with the dE VCC 
1947. The two systems compare quite well for the observed 
part. Figur^3]bottom panel shows the surface brightness profile 
of the initial model which can be fitted to a Sersic profile with 
index n — 1 .07, similar to dwarf galaxies in the Coma cluster 
(Aguerri et al. 2005). 

The scaling used (M - r 2 ) implies a constant surface den- 
sity for all our dwarf models. This is at odds with the relation 
between surface brightness and magnitude observed for dwarf 
galaxies (Ferguson & Binggeli 1994). However this relation 
and the one found by Aguerri et al. (2005; see their Fig. 8) show 
a large scatter that allow us to consider a constant central sur- 
face brightness for several magnitudes as a first approximation. 
Our initial conditions cover 4 magnitudes in absolute magni- 
tude, so for our smaller model CeP\ we should either consider 
a different M/L to comply with such relation or build the ini- 
tial model using a different scaling. For consistency reasons we 
choose the first option, a discussion on the effects of the second 
one is given below. 

In order to test the effects of a dark matter halo around the 
dwarf, we also run models with a dwarf consisting of a King 
model inside a halo. The halo is an Evans model (Kuijken & 
Dubinski 1994); the global central potential is -4.6, and the 
halo asymptotic circular velocity is 2. These parameters give a 
total mass-to-luminous mass ratio of M/L m 8, which is on the 
low regime of values for the dwarf galaxies in the local group 
(Mateo 1998). 

The more massive galaxy is simulated as an isotropic, non- 
rotating Jaffe (1983) model. The theoretical half mass radius 
O and the total mass M\ were set to unity, so that we have a 
mass ratio for the encounter of 10:1. The initial system was 
let to evolve in isolation for 10 crossing times. After a mild re- 
organization, the system reaches equilibrium. Because a cut-off 
radius was imposed at r - 10 xrj, the true half-mass radius r\/2 
is 0.83. 

Non-dimensional units are adopted throughout, with the 
constant of gravity G=l. The theoretical half mass radius of 
the Jaffe model o and the total mass of the elliptical galaxy are 
also set to unity. A set of physical units to compare our models 
with real galaxies would be the following: 

[M] = M\ = 4 x 10 11 M Q , (5) 

[L] = n = 10 kpc, (6) 

[T] = 2.4 x 10 7 yr. (7) 
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Fig. 3. Top: rotational support, expressed as V(R)/cr(R) vs. 
R/Re of the initial model (solid line) compared with the dE 
VCC 1947 (diamonds). Bottom: Surface brightness profile of 
the initial model. 



With these, the velocity unit is: 
[v] =414km/s. 



(8) 



With these units the initial maximum rotation velocity of the 
dwarf galaxy is « 60 km/s. 

Initial conditions for the fly-by experiments are listed in 
Table [0 Column (1) gives the model code name. Column (2) 
lists the mass ratio. Column (3) gives the spin orientation of the 
smaller galaxy. Columns (4), (5), (6) and (7) give the initial sep- 
aration, orbital eccentricity, pericenter distance and pericenter 
velocity, respectively, of the orbit the galaxies would follow if 
they were point masses on Keplerian orbits. The initial separa- 
tion was always equal to 20 length units, or twice the radius en- 
closing 99% of the mass of the more massive galaxy. Inspection 
of Table[2shows that most of our models explore the regime of 
parabolic or hyperbolic orbits with small pericenter distances 

^"peri ■ 

Spin orientation is coded by the first upper-case letters in 
the model names. C denotes a counterspinning dwarf. Given 
that we are searching for counterrotating signatures, most of 
our models are C models. However, we did run one model with 
dwarf rotation aligned with the orbit (denoted with first letter 
R) and one model with partially misaligned dwarf spin (first 
letters CR). The lower-case letter in the model name codes r pe[ j, 
which varies from r per i= 1/4 of the half mass radius of the dwarf 
galaxy (models a) up to a r per j equal to the radius including 90% 
of the mass of the elliptical (models /). For each r per j we run 
hyperbolic (H, e = 1.07) and parabolic (P) orbits. For model a, 



the P orbit results in a fast merger, and is not described further 
here. For models a and b, we additionally run hyperbolic orbits 
with e = 1.003 and e = 1.043, which have a longer interaction 
time (labeled with HH). In general the time it takes to the dwarf 
to cross the inner regions of the giant is less than 1 half-mass 
crossing time of the giant. 

Models listed in the lower part of Tableware variations on 
the models listed above, in which we selectively change spe- 
cific model or orbital parameters. These models are intended to 
provide clues on the effects that a given parameter may have on 
the generation of core-envelope counterspinning. Model CePc 
is a re-edition of model CeP using five times more particles; 
important differences between CeP and CePc would point to 
resolution issues in our results. Model CeP e again re-edits 
model CeP using a smaller value of the softening parameter 
in the gravitational calculation (see below). The effect of a dark 
matter halo surrounding the dwarf galaxy was investigated with 
model CePii which again has identical orbital conditions to 
CeP. In models ReP and CReP, the spin arrangement of model 
CeP is modified to test the effect of a partial alignment of the 
secondary spin with the orbit. Three models were run with dif- 
ferent mass ratio and similar orbital conditions to CeP (mass 
ratio 5:1, 20:1 and 50:1 for models CeP2, CeP3 and CeP4 re- 
spectively). 

We also run three more models to study the influence of 
the initial rotation curve of the dwarf galaxy in the final results. 
The initial dwarf galaxy for these models was the final state of 
the dwarf from model ReP. We have placed it on an parabolic 
orbit with anti-parallel spin and r pei i = 0.83, 0.2 and 0. 1 . These 
models are named CePu, CcPr, and CbPu, respectively. 

We used the tree-code GADGET 1.1 (Springel et al. 2001) 
to run our experiments on a Beowulf cluster making use of 
16 CPU's. We used 100000 particles for the King model and 
102400 particles for the Jaffe model. Individual softening for 
each component was used with e ~ 1/5 of the half mass radius 
system (e = 0.04 for the dwarf and e = 0.1 for the elliptical; 
the check-up run Cef e uses values of s = 0.015 and e = 0.045 
respectively). The tolerance parameter was set to 9 — 0.8 and 
quadrupole terms were included in the force calculation. A typ- 
ical run with 2 x 10 5 particles took of the order of 1.8 x 10 4 
seconds for 7 x 10 3 steps. Conservation of energy was good, 
with errors well below 0.1%. 

All models were let to evolve well after the first pass 
through the pericenter, until the distance between the cen- 
troids of the two systems was roughly the initial separation. 
Dynamical friction prevented models CaHH and CbP from 
even approaching the initial separation; in these cases, a dis- 
tance close to apocenter after the first pass through the peri- 
center was taken for analysis. Rotation curves for the dwarfs 
were derived by observing the galaxy models from a point of 
view perpendicular to the orbital angular momentum and also 
perpendicular to any tidal feature that might be present in the 
final dwarf galaxy. Variations on the amplitude of the rotation 
curve are to be expected when viewed from different points of 
view. The rotation curves were obtained by placing a slit along 
its projected major axis with a length of at least 2 length units 
(=20 kpc, 20 x ri/2 of the dwarf galaxy) and a width of 0.3. 
They give mass-weighted line-of-sight velocities. 



6 



Gonzalez-Garcfa, Aguerri and Balcells: Harassment Origin for KDC in dE's? 






Fig. 4. Rotation velocity curves for the different models, (a) 
Models with hyperbolic energies and pericenter distance less 
than R e of the dwarf galaxy, (b) Models with hyperbolic ener- 
gies and pericenter distance larger than R e of the dwarf galaxy, 
(c) Models with parabolic energies 



Fig. 5. Rotation velocity curves for check up models, (a) 
Models starting from different initial conditions, (b) Models 
with different initial mass ratio, (c) Model with 5 times more 
particles, model with a halo and model with a different soften- 
ing. 



4. Results 

Rotation curves for the models from Table 1 are shown in 
Figures |4] and [5] Models are grouped in the following way: 
Fig. ^fr,b: hyperbolic encounters, sorted by increasing r per i. 
Fig.^p: parabolic encounters. Fig. [5^: variations of the dwarf 
orientation and rotation curve shape. Fig. [5j?: variations of the 
mass ratio. changes in the number of particles, and on 

the presence of a halo in the dwarf. In each of the panels, the 
rotation curve for the initial dwarf model is shown for compar- 
ison (solid line). 

Figures |4] and [5] show that the rotation amplitudes of most 
models are lower than in the initial model: as expected, the fly- 
by interactions do produce a transfer of orbital angular momen- 
tum to the dwarf. For retrograde experiments with small r per j, a 



counter-rotation (CR) signature appears. However, a quick in- 
spection reveals that the radius of rotation curve reversal R cl is 
significantly larger in the models (~ 2-3 R e 2) than in the ob- 
served galaxies. R CI decreases monotonically as r per j decreases. 
In order to have R CI ss R e2 , the encounter has to be almost head- 
on and mildly energetic (runs CbP and CaHH). 

Figure 0? focuses on the interesting regime of hyperbolic, 
nearly-head on collisions. Two values of orbit eccentricity and 
two values of r per j are shown. For the two most energetic mod- 
els (CaH and CbH) we find R a /R e 2 > 3, showing that nearly- 
head on orbits do not ensure that a CR signature appears at 
radial distances similar to those of real dwarfs. Angular mo- 
mentum transfer to the dwarf has low efficiency because the in- 
teraction time in these fast encounters is short. A more realistic 
CR radius (R/R e 2 ~ 1.5) is obtained in the models with lower 
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eccentricity {CaHH and CbHH). Hence, the parameter space 
or orbits leading to CR signatures in dwarfs is constrained both 
in r per ; and in orbital energy: fast orbits yield too short an inter- 
action time, while orbits approaching parabolic conditions lead 
to a fast merger. 

We have checked for the stability of the counterrotation by 
letting the remnant of model CaHH to run in isolation for 50 
time units (equivalent to 12 Gyrs). The counterrotation is stable 
throughout this time. 

5. Dependence of the results on the encounter 
parameters 

The results from the previous section provide measurements 
on how impact parameter and encounter speed determine the 
onset of envelope CR in fly-by interactions. Despite the limited 
number of models we have run, the results give clues that only 
very strong interactions are capable of imprinting CR at radii 
similar to those observed in the galaxies FS373 and FS76 (§0. 
We now explore to what degree the results change when other 
interaction characteristics vary. It is not our purpose to do a full 
exploration of parameter space; we just replicate some of the 
models discussed in §@] varying one single parameter in turn. 
We have changed: spin orientation, initial rotation curve of the 
dwarf galaxy, mass-ratio, presence of a dark matter halo in the 
dwarf galaxy, particle number and softening. Initial conditions 
for these models are listed in Tabled We discuss these models 
in the following subsections. 

5.1. The spin orientation 

Three of our models have differing spin orientations, while 
keeping all other parameters unchanged: CeP (dwarf spin an- 
tiparallel to the orbit), ReP (spin aligned with the orbit) and 
CReP (spin at 135° to the orbit). When the dwarf spin is aligned 
with the orbit (ReP, Fig: |5^), the external parts of the dwarf 
gain angular momentum and rotate faster than initially. The fi- 
nal rotation curve is almost flat. No counter-rotating features 
appear in this model, as expected from theoretical expectations. 

The final system for model CReP is quite similar to the 
model CeP, with R a /R e 2 ~ 3 (compare with Figs.@J; and|5^). 
This result suggests that the onset of CR does not depend on 
the dwarf spin being precisely antiparallel to the orbit. 

5.2. Shape of the dwarf rotation curve 

As a result of the method employed to impart rotation (§0, our 
initial dwarf galaxy model has an outwardly declining rotation 
curve. This is in contrast to rotation curves of dwarf galaxies, 
which typically increase with radius (e.g., Swaters et al. 2003; 
van Zee et al. 2004). As a simple check on the effects of this dif- 
ference on our results, we have run experiments using as initial 
conditions the dwarf galaxy after the interaction in model ReP. 
As shown in the previous subsection, the rotation curve after 
fly-by ReP is shallower in the inner parts than our canonical ini- 
tial model, and is nearly flat in the outer parts (see Fig.[5k). We 
have run 3 models with different r pe i-i : CbPu, CcPr and CePu. 



The results are shown in Figure^. Again, we find counter- 
rotation, but R cr of these new models is larger than for those of 
models CbP, CeP and CeP. This is due to the larger rotation 
in the outer-most region of the new initial model. These exper- 
iments also show the possible cumulative effect of harassment. 
After two encounters the CR appears at larger radii than for a 
single encounter. Thus, in principle, repeated encounters would 
make even harder to find CR. 

5.3. Models with different mass ratio 

Models CeP2, CeP3 and CeP4 use the same initial orbital pa- 
rameters as model CeP, but we modify the mass ratio, which 
is set to 5:1, 20:1 and 50:1, for models CePl, CeP3 and CeP4, 
respectively (Table^. The final rotation curves from the three 
models are shown in Fig.[5J). R cr /R e 2 is quite similar for the two 
high-mass satellite cases, but is significantly larger for the less 
massive dwarfs (R CI /R e 2 ~ 4). Clearly, smaller mass ratios fa- 
vor the appearance of counterrotation at smaller radii. We note 
that R cl /R e 2 increases for larger mass ratios because the effec- 
tive radius of the smaller system decreases faster than the R a 
does. This is due to the difference in densities introduced by 
the physical scaling used in our models. Smaller systems are 
denser in the inner parts and this prevents against the forma- 
tion of the counter-rotation in the inner parts. Extrapolation of 
this trend suggests that fly-by interactions might yield smaller 
values of R a /R e 2 when applied to the less extreme mass ratios 
typical of encounters between dwarf galaxies. 

5.4. Presence of dark matter halo in the dwarf galaxy 

We have run a model with the same initial orbital parameters as 
CeP but including a dark matter halo around the dwarf galaxy 
(model CePn). Results for this model and model CeP are com- 
pared in Fig.|5J;. We find that model CePn shows no counter- 
rotation over the explored radial range. Clearly, the dark matter 
halo of the dwarf galaxy absorbs a large fraction of the orbital 
angular momentum leaving no significant counter-rotating sig- 
nature on the final system. The halo-embedded dwarf CePn 
has a shallower inner rotation curve after the interaction, point- 
ing at an enhanced ability to retain orbital angular momentum 
thanks to the deeper potential well provided by the dark halo. 

5.5. Models with different number of particles and 
softening 

In order to check for resolution problems affecting the results, 
we have run a model which has the same initial conditions as 
CeP but with 5 times more particles (model CePo). The results 
are shown in Fig.|5J;. We find that the rotation curves of models 
CeP and CeP c are identical, suggesting that particle number 
does not critically influence the results presented here. 

We have also run model CeP with a smaller softening in 
order to see the effect of having a more accurate force resolu- 
tion. The results are shown in Fig.|3J;. We find that the rotation 
curves of models CeP and CeP e are similar but not identical. 
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Fig. 6. Angular momentum change versus pericenter distance. 
The theoretical expectation from the pure impulse approxima- 
tion for our runs are given as open circles (see eq. |9} and the 
full approximation is given as plus signs (see eq.H0>. The val- 
ues of the angular momentum change for our models at radii 
including the 50, 70, 90 and 100% of the mass of the model are 
given as diamonds, triangles, squares and crosses respectively. 
The shadowed area shows the region where the pure tidal ap- 
proximation applies. The thick and light arrows indicate that 
AJ tends to zero for head-on encounters (r per i= 0). 

The counterrotation in the new model appears for outer 
radii than in the fiducial model. 

6. Discussion 

6.1. Numerical models vs. impulsive approximation 

Comparison of the rotation curves from our Af-body fly-by in- 
teractions (Figs.|4]and|5} with the theoretical expectations (§|2] 
and Fig.Q shows that the fly-by interactions are far less effec- 
tive in reversing the velocity sign than expected from the pure 
impulse approximation. The full treatment must be applied in 
order to account for the observed behaviour. 

A first clue to this behaviour is given by the radial distri- 
bution of A(J 2 ), the angular momentum deposited in the sec- 
ondary. In Figure[6]we plot A(/2) integrated out to radii includ- 
ing 50%, 70%, 90% and 100% of the secondary mass, against 
the pericenter distance. The figure indicates that, for all the 
models, A{J 2 ) in the outer 10% of the mass is over 10 times 
larger than that in the inner 90% of the secondary mass. As 
suggested by Figure |2 the outer parts of the satellite become 
unbound after the interaction. Therefore, over 90% of the trans- 
ferred angular momentum escapes the dwarf in tidal tails. The 
effects of the interaction on the dwarf rotation curve must be 
weaker than predicted by eqn.^ 

A second clue to the discrepancy is provided by comparing 
the true angular momentum transfer to that predicted by the 
pure impulse approximation in the tidal limit. Using eqn. ^ 
AJ 2 (R) ~ M 2 RAV 2 (R), and I 22 (R) ~ M 2 R 2 /3, we get: 

AJ 2 (R) ~ 2^-I 22 (R), (9) 



where I 22 is the smallest eigenvalue of the inertia tensor. A(J 2 ) 
computed with eqn. [9] for the entire secondary is plotted with 
circles in Figure [5] We find that eqn. [5] accurately predicts the 
total angular momentum transfer (crosses in Fig.|6jl for exper- 
iments with intermediate pericenter distances, but it overesti- 
mates A(J 2 ) at small r per i, and it underestimates A(J 2 ) at large 
''pen- The overestimate of A(J 2 ) at small r per i aggravates the dis- 
crepancy between the prediction of the impulse approximation 
and the fly-by experiments. As we diminish r pe ri, A(J 2 ) tends to 
zero, as we approach head-on conditions (indicated with arrows 
in Figure^. At the other end, log(r pe ,i/^ e 2) ^ 1, the discrep- 
ancy between measured and predicted A(/2) may be due to the 
fact that the interaction time is comparable to the dynamical 
time of the secondary, in contradiction with the assumptions of 
the impulse approximation. 

Figure [6] supports the findings of Aguilar & White (1985), 
who show that the tidal approximation gives reasonable results 
for impact parameters of order 5 x r/,, with r/, the median radius 
of the small system. This range of impact parameters, shaded 
in figure^ encompasses the models where eqn. [9] reproduces 
the measured A(J 2 ). 

The discrepancies are less extreme if we take into account 
the full treatment of the impulse approximation, and use eqn. [2] 
Then we get: 

AJ 2 (R) = 2^/ 22 (tfW-^)(l - ojt)- 2 - 5 (10) 

This is given in Figure [6] as plus signs. The agreement is quite 
good even for the inner parts. For large pericenter distances the 
approximation does not apply due to the long duration of the 
tidal impulse and we are still not able to account for the angular 
momentum carried away by the tidal tails. 

6.2. Comparison with observations 

As described in §[T] de Rijcke et al. (2004) have found KDCs 
in two dE elliptical galaxies belonging to the NGC 5044 and 
NGC 3258 groups. These galaxies show CR features at R a /R e ~ 
1 — 1.5, which de Rijcke et al. explain as an effect of harassment 
interactions. Other observational studies of dE rotation curves 
rarely reach the depths required to map rotation to ~R e , hence 
determining the frequency of CR in dEs is subject to strong ob- 
servational biases. The studies of Virgo dwarfs (Pedraz et al. 
2002; van Zee et al. 2004) do show a fraction of objects with 
kinematic substructure, such as asymmetric or otherwise com- 
plex rotation curves. In particular, 1 out of 6 dEs from Pedraz 
et al. (2002) shows hints of envelope counterrotation. 

Figure summarizes our results showing the counterrota- 
tion radii vs. the pericenter distance in units of the secondary 
effective radius. The values for the CR in the two dE from de 
Rijcke et al. (2004) are given as a comparison. Interactions such 
as our models a and b may explain the CR features observed in 
the two dwarf galaxies studied by de Rijcke et al. (2004). 

However, our models tend to produce CR signatures at 
larger radii than is observed. The required small impact param- 
eters are a strong restriction for the types of collisions that may 
lead to counterrotation. In the next subsection we estimate the 
probability of such collisions in a cluster environment. 
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bital configurations. 

6.3. The harassment scenario for dwarf CR 

To determine whether harassment is a likely mechanism for 
the formation of CR, we estimate here the probability of small 
impact parameter collisions. Rather than a group, we choose 
a cluster environment since it contains large numbers of giant 
and dwarf galaxies, which allows for a statistical calculation. 
We use a very simplified picture of the dynamical state of a 
cluster, which we describe with a single, representative num- 
ber density of large galaxies, and we assume that velocities 
are uncorrected. The cross-section E CR for collisions leading 
to core-envelope CR is approximately 

£ CR ~ 0.15tt£ 2 (11) 

where b » 5 kpc is a typical impact parameter for experiments 
Ca and Cb, and the numerical factor is the fraction of encoun- 
ters with spins within 45° from anti-parallel. The mean-free 
path A CR between collisions leading to CR is 



where n is the number density of bright galaxies in the cluster. 
For Virgo, taking n » 50 galMpc~ 3 (Binggeli et al. 1987), we 
infer A CR w 1700 Mpc. In a time interval At, a dwarf in a cluster 
with one-dimensional velocity dispersion <x 1D will travel a dis- 
tance L = V(3) tr 1D A t, hence the probability of a CR collision 
may be expressed as 

Pgr = L/A CR = V3 cr ID At n 2 CR , (13) 

then, for N D dwarfs, the number of them prone to have a CR 
collision is 

N C r =N D V3o- 1D Af«2 CR , (14) 

which, for values typical of the Virgo cluster (Binggeli et 
al. 1985), yields 



N C R = 




In other words, we may expect to find 8 dwarfs showing CR 
in a cluster like Virgo. This means that, at face value, about 1% 
of dwarfs in a Virgo-type cluster may experience a CR collision 
in 10 Gyr. Our estimate is admittedly crude. Some of the col- 
lisions will be too slow, or too fast, for an efficient transfer of 
angular momentum; dwarfs embedded in their own dark mat- 
ter halos hardly respond to the interaction; internal velocities in 
the cluster are not uncorrelated; and, dwarfs in elongated orbits 
spend most of their lifetimes in the low-density outskirts of the 
cluster. On the other hand, the cross-section for CR is bigger in 
interactions with a smaller mass ratio than modeled here if, as 
suggested in § 15.31 lower mass ratios lead to smaller R CI /R e 2- 

Finally, more complex prescriptions for the structure and 
internal dynamics of the dwarf may lead to different results 
on R cl . For instance, a dwarf consisting of a disk embedded 
in a non-rotating envelope requires less angular momentum 
exchange to generate a CR signature, as shown by the fly-by 
models explored by Hau & Thomson (1994). The scaling used 
for our initial conditions (Fish 1966) implies a constant surface 
density. As said in section^] this could be consistent for a range 
of magnitudes for the dwarf galaxies studied by Aguerri et al. 
(2005). A different scaling may reproduce the general trend ob- 
served where fainter galaxies have a lower surface brightness. 
This would imply that the volume density would not increase 
with decreasing mass as fast as in our scaling, changing some- 
what the CR radius for the smaller systems. 

Despite these caveats, our conclusions hold and in a cluster 
like Virgo ~ 1 % of dwarf galaxies may present counter-rotation 
formed by harassment. 
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